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ABSTRACT 



^Or- 

v//^ '-5 ,> 

A control system was designed to attenuate vertical ^ 
accelerations for the XR"3 captured air bubble type surface 
effect ship using linear regulator techniques applied to the 
simplified nonlinear equations of motion. A pressure lift- 
only model was used to represent the craft vertical heave 
motion and was linearized around the steady state operating 
point. Model validation was obtained through analysis of 
the frequency spectrum. State variable feedback was used to 
determine a set of optimal control gains that would reduce 
the magnitude of the heave acceleration during operation 
under simulated sea input conditions. 
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I , INTRODUCTION 



An important design consideration for the high speed 
captured air bubble type surface effect ship is the vertical 
acceleration associated with operation in a random ocean 
environment [Ref, 1], The vertical motion, heave, is caused 
by waves changing the volume of the plenum producing pres- 
sure variations that result in vertical craft accelerations. 

Studies of conventional hull ships show that the strongest 
vertical plane accelerations are experienced at encounter 
frequencies close to the ship’s natural frequency. Reference 
2 states that "The most severe and objectionable motions (wet 
decks, slamming, high accelerations) result from the responses 
to the components that are near synchronism with the ship's 
natural pitching and heaving periods" . Encounter frequency 
is a function of wave frequency, craft speed and direction. 
Since most of the weight of the surface effect ship is sup- 
ported by a bubble of air, large accelerations can occur that 
affect habitability and sea keeping ability while operating 
in high seas . 

Methods of attenuating vertical motion have been the 
subject of past studies. A conventional design by Rohr Marine 
Inc. [Ref. 3] utilized accelerometer output to monitor verti- 
cal motion variations. These signals are passed through 
compensation circuits, filtered and applied to actuators 
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which control air flow rates into the plenum. State variable 
techniques are applied to describe the system and control 
functions. The physical constraints of the actuators are 
considered with regat'd to their inherent nonlinearities and 
limitations , The sea input is modeled using a wave height 
power spectral density procedure. 

Through the use of a towing tank model and data gathered 
from the 10 OB test craft, an investigation was conducted to 
measure the characteristic response of the SES while operating 
in an actual sea state environment [Ref. 4], The model was 
subjected to simulated sea state inputs of various severity in 
an effort to establish the linear characteristics of the 
vertical motion and determine the damped natural frequency of 
response. A series of curves showing response amplitude for 
both C.G. accelerations and pitch angle versus encounter fre- 
quency were developed for a range of speed and sea state 
operation . 

In Ref. 1, the frequency response approach was used in an 
analysis which examined the effects of changing plenum air 
flow rates on center of gravity accelerations using the six 
degree of freedom simulation program of the XR-3 [Ref. 5]. 

This study developed a linear simplified model in close agree- 
ment with the 6 D.0,F. model and in addition used signifi- 
cantly less computer time, 

In Ref, 6, Boggio conducted a study on the effect of the 
use of a flexible membrane in the plenum chamber itself. 

Using the loads and motions program for the XR-3, the study 
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reduced computation time through the inclusion of pressure 
rate and volume rate equations, Membrane insertion served to 
effectively reduce the magnitude of pressure variations in 
the plenum thereby attenuating vertical motions. 

Grant [Ref, 7] used simple models for the control devices 
applied to the loads and motions program of the lOOB SES 
craft. From monitored plenum pressure variations and heave 
accelerations, control signals were developed by filter de- 
sign methods , to vary the fan speed and vent louver openings 
in order to reduce the vertical accelerations and pitch motion. 

This thesis applies the linear regulator design technique 
to determine a set of optimal state feedback gains to control 
the air flow rate . The equations of motion used in the design 
were derived from a simplified pressure lift only model lin- 
earized around the steady state operating point. Design re- 
sults were evaluated using both the linear and nonlinear 
models . 
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II. PRESENTATION OP THE SYSTEM 



Two major forces which support the weight of the surface 
effect ship are the buoyancy and plenum pressure forces. 

The buoyancy forces are exerted on the craft by dis- 

placement of water by the underwater volume of the rigid 
sidewalls. The plenum pressure lift force (Fp) is exerted 
on the craft by the captured cushion of air. Flexible bow 
and stern seals along with the rigid sidewalls and lift fans 
form the captured air bubble. Relative contributions to 
craft vertical plane motion by Fj^ and Fp depend upon the 
encounter frequency • At low frequencies, Fj^ responds 

to the rise and fall of the sea. At high , Fj^ attenuates 
because incremental volume changes of the immersed hull are 
reduced. The Fp force responds inversly to changes in plenum 
volume (V„) according to the adiabatic law of pressure 

D 

variations within the plenum chamber [Ref. 1]. Volumetric 
changes in the plenum depend upon the sea state and encounter 
frequency. As is shown in Ref. 1, the effect of the volume 
changes produces an increase in the Fp force as the encounter 
frequency is increased. 

In the system developed in this study, the buoyancy terms 
have been omitted in an effort to simplify the system further. 
This simplified model is referred to as the ’’pressure lift 
only" system. Transfer functions of acceleration, volumetric 
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plenum air mass to input waveheight and system state equa- 
tions were developed from the linearized form of the pres- 
sure lift only system. 

The basic assumptions which form the basis for the pres- 
sure lift only model are presented below: 

1. The center of pressure is directly under the 

center of gravity (C.G,), 

2. The sidewalls are of uniform cross section and are 
symmetrical about the C.G. 

3. The effect of pitch moments on the vertical plane 
forces are neglected. 

4. The only variable lift forces are those forces 
caused by changes in plenum pressure (Pp)* 

5. The craft is at a constant speed. 

6. The air leakage rate from the rear seal is constant. 

7. The craft weight is supported by two component 
forces in the following approximate proportions ; 

Sidewalls: 20% (fixed force) 

Captured Air Bubble: 80% (variable force) 

8 . Encounter frequencies of interest lie within the 
range of 17<w^<100 radians/sec. 



A. EQUATIONS OF MOTION 

The development of the nonlinear simplified model for 
vertica.1 motion is presented in Ref, 1 and the block diagram 
appears in Appendix A. The development of the pressure lift 
only nonlinear model follows in a similar manner. The equa- 
tions of motion differ in that the force is considered 
constant . 
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The application of the given assumptions to the Newtonian 
Laws for vertical motion of a rigid body system leads to the 
following equations; 

1) Orifice Leakage Rate 



q , = C A, 
^out n 1 



Vf 



£11 

sec 



(2.A.1) 



2) Fan Map Input Rate 



q.„ = n[Qi<o) - 



ft 

sec 



(2. A. 2) 



3) Absolute Plenum Pressure 



M Y 

P, = P ( rf-^ ) — 7 

b a V, p. 2 



lbs 



b^a 



ft 



(2. A. 3) 



4) Plenum Air Flow Rate 



lbs 

b ^ '"a "^in ~^out' sec 






(2. A. 4) 



5) Heave Acceleration 



•• T7 -F-h 2 

7 _ /W _P _!•) It 

^"^M"M~M'^ sec 



(2. A. 5) 



6) Heave Velocity 



Z = 






dt 



ft 

sec 



(2. A. 6) 

% 
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7) 



Heave Displacement 



Z = y* Z dt 



ft 



(2,/i 



8) Plenum Pressure Lift Force 



"p - 



lbs 



(2. A 



9) Buoyancy Lift Force 



Ft = [2A l,p ] 
1 s 



lbs 



(2 .A 



where 1^ is steady state draft, 



10) Plenum Gage Pressure 



P, = (P- - P ) 

b b a 



lbs 

ft^ 



(2. A 



11) Draft 



1 , = Z + Z 
d s 



ft 



(2. A 



The system parameters and constants are listed as follows: 



1) Adiabatic Process Coefficient 



Y = 1.4 



, 8 ) 



.9) 



. 10 ) 



. 11 ) 
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2) Leakage Area 



= ,438 ft^ 



3) Leakage Orifice Coefficient 



C = .90 
n 



4) Air Density 



p = .002378 



5) Atmospheric Pressure 



P = 2116 
a 



6) Plenum Area 



A, = 200 

D 



ft 



2 



7) Empty Plenum Volume 

V = 383 ft^ 

n 

8) Craft Weight 



W = 6007 lbs 
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9) Area of Each Inunersed Sidewall 
A = 75/4 = 18.75 ft^ 

10) Draft Initial Condition 

L^(o) = ,5 ft 

11) Density of Water 
p = 1.99 

12) Gravity Acceleration Constant 

G = 32.174 ft/sec^ 

13) Vertical Center of Gravity Location Above Keel 
Line 

Z = 2.5 ft 

s 

14) Number of Fans 
n = 5 

15) Steady State Fan Output 

Q . (o) = 35 ft^/sec 
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16) Fan Map Slope 



K =1.0 

q 
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0) 



3: 

v_^ 

f“ 



fH na 0) 
0) *H W 

I— ! \ 

O fU 
> fH 



0) 

D 

bO 

•H 
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PRESSURE LIFT ONLY SCHEMATIC DIAGRAM 



B. LINEAR MODEL 



Prior to defining the system sta.tes, the non linear 
equations were linearized around their steady state operating 
conditions by application of a Taylor series expansion. The 
details of the linearization process is shown in Appendix A. 

The resulting simplified pressure only signal flow graph 
for the linearized system is shown in Figure 2. 

1 , Transfer Function Development 



The system graphical determinant is; 



A = 1 - (L^ + L2> + L^ L^ 



(2.B.1) 



A = 1 - (L^ + 



(2 .B. 2) 



where 



L^ = s”^s~^(-kg) (-k2) (k^) 




(2.B.3) 



L_ = s~^(k_) (kj(-k-) 
2 7 0 6 




(2.B.4) 



The transmission gain from the control input to the accel- 
eration output for calm water operation is ; 
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Figure 



(2.B.5) 









a 



u 



where 



^1 ■ 



-1 



(2 .B. 6) 



and 



A 




(2.B.7) 



and 



a 



w 



0 



(2 .B. 8) 



Substituting equations (2.B,2), (2.B.6), (2.B.7) into 
equation (2.B.5) yields; 



3“ 

z 



-k3k^k7kj^^,s'^ 8u 

l-C-s’^k^kjk^^kg ) + C -s'^k^kgky)] 



(2 .B. 9) 



Reducing further; 



“^ 3 ^ 4 ^ 7^10 



u 



'^z = -T 

s +]c 2 J<:gl <^7 s+k 2 k 3 k^^kg 



(2.B.10) 



The transmission gain from the sinusoidal sea generator input 
to the a,cceleration output with no control signal is; 






^ 2^2 



"W 



(2.B.11) 
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where : 






(2,B.12) 



and 



= 1 (2.B.13) 

and 

3^ = 0 (2.B.14) 



Appropriate substitution into equation (2.B,11) 



yields 



1-C (s~^k2k2k^kg )+( -s'^k^kgky ) ] 



Reducing further, the transfer function is 



a 

9W 



2 

-k^k 2 kgk^^ s 
s +k 2 kgk 2 S+k 2 kgk^^kg 



(2.B.15) 



For the condition that control input is applied, the 

acceleration output has the following transform; 

-k„k,,k„k, 3 

3 4 7 10 u 



s +kgkgk 2 S+k 2 kgk^^kg 



^ 1 ^ 2 ^ 3 ^ 4 ^ 



9W 



s +kgkgk 2 S+k 2 kgk^kg 



(2 .B.16) 
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Following the same procedure for the applied control 
input, the volumetric plenum air mass output has the 
following transfer function; 



43 = 1... 



^3 = '"10® 



"1 



A 2 = 1 



(2,B,17) 



(2.B.18) 



(2.B.19) 



-k^k^k^kg 



9V 



m 



s^+kakakTS+kjkjkjjkg 



U 



Equation (2,B,15) is of the form: 



(2.B.20) 

(2.B.21) 



OUTPUT 

INPUT 



2 

-as 



s^+bs+c 



(2 .B. 22) 



and with no control action, becomes*, 



where 



y(s) _ as' 






(2.B.23) 



y(s) = 9Z 



(2 ,B. 24) 



and 



, b 
s^+bs+c 



(2 .B, 25) 
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where the coefficients consist of the following parameters j 



a -= 






b = 



kskgk^ 






d =. 



-k^k^kgkg 



(2.B, 26) 
(2.B.27) 
(2.B.28) 
(2.B.29) 



The values of the coefficients, k^, for the XR-3 craft are 
listed in Table I and were obtained from the results of the 
linearization process of Appendix A, 

Special attention should be given to the coefficient 
that contains k^ which is a sine function of the form; 

2W V w.L 

V = Sin (2.B.30) 

1 to. 2VW 



k^ may be simplified by an approximation detailed in Appendix 
A of Ref. 8. The simplification is based upon the fact that 
for ahead seas, the incident frequency, to^ is-, 




1 “• y 1 1 4 ( ) 0) 

-2(-#;) 



(2 ,B. 31) 



For convenience in this study, the value of the ships velocity, 

V is set equal to the gravitational constant, G. This value 
s 

of V , when converted to knots, is 19.3 knots which is within 
s ’ 



23 



TABLE I 



SYSTEM COEFFICIENT (k ) VALUES 

n 



FORMULA 
2WpV„ W L 



03 . 






w 



1 

V, (0) 

D 



XR^3 VALUE 
2267,782 



sinC 



62 .8319 



/62.8319s 

W 



) 



3 ,5336x10"^ 



YPb(O) 



2996 . 0 




nk 



q 



C A, / P 

-g-i / — 5~ 

a 2P7(0) 

D 



1,0711 
6 , 14422 



1 

V (0) 



m 



3 .5052x10 




200.0 



CONSTAiNT 



1.0 



the operating range of the XR-3, Continuing; 






i 




C2.B.32) 



It follows that 



u 



2 

i 




(2.B. 33) 



This approximation is valid for all tj^>2,5 radians/sec . From 
the definition of the incident wave , 



u)i 



2'ITg 

Hi 



(2.B.34) 



becomes 



■“iS 



Sin" 2V 



= A. 



hS 



w.L 



2V 



tu 



(2.B.35) 



Where the horizontal area of the plenum chamber is 



Ab = WpLp 



(2.B, 36) 



For a speed of 19.3 knots and since 



L 



W 




C2.B.37) 
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simplifies toj 



TfL, 



Sin(y.^) 

\ ""Ju ' 



For the XR'->3 , 



A, = 2 00 ft' 

D 



L = 20 ft 
P 



Then 



a = = 2267,782 



Sin(i2;iB9) 



.62.8319. 
^ "L.: 

W 



From Equation (2.B.38), it can be seen that a null 

k, will occur for nL„ = L and a peak value for k, 
1 W p ^ 1 

for nL^^ = Lp/2 where n is an integer. 

2 , State Equation Development 

The state equations are of the form 

X = AX + B3u + G3W 



where the state variable vector is chosen as 

X = O , *3’- , 3V 
— z ’ z ’ m 



and the output variable 



y=3 = CX + F3W 

z f, - 



(2.B.38) 



(2.B.39) 

value for 
will occur 



(2.B.40) 

(2.B.41) 



(2 ,B.42) 
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0 



ir 






EW 



I 



The states are defined as 



XI = 


9Z 


II 

CSl 

) 

tsi 


X2 = 


9Z 


= (Z-Z 


X3 = 


9V 


= (V 



m m 



ss 



)> 



ss 



), 



-v, 



mss 



Incremental Heave (Draft) 

Incremental Heave Rate 

) , Incremental Air Mass 
Volume 



where V 

m 




P 



a 



with control defined as a scaler, rate of incremental 
of air mass volume, 



3 



u 



9V 



m 



c 



and where the sea input disturbance is 9W. 

From the signal flow graph of Fig. 2, 



9V = (9 k^^=k-9P, )s 
m u 10 6 b 

X^. = ~ s"^9z 

• , —1 

X, = 9^ = s 9z 

1 z 



X 



2 



3Pb^ 



(2.B,43) 

C2.B.44) 

(2.B.45) 

(2.B.46) 

change 

(2 .B.47) 

(2.B.48) 

(2.B. 49) 
(2 ,B. 50) 

(2 ,B, 51) 
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where 



Rearranging 


C2.B, 52) 


8P, = av k,k,,+k, k„k^ 3W+k„k_k„ 9Z 
b mo / 123 238 

and grouping in order of states 


(2.B, 53) 


3P, = k-k.k„9Z+k»k„ 3V +k^k.k,9W 

b 238 37ml23 

The output term, y is 


(2.B.54) 


II 

tsi • 


( 2,B. 55) 



Substituting equation (2,B,5) into equation (2.B.55) yields, 

y=-k.k-k,,kQ8Z-k.k„k„aV -k,k.k,k„9W (2.B.56) 

zc54o J4 / m lzo4 

Finally, the incremental change in the volumetric air mass 
is written as 



X3 = aV =( auk, n-k. 3P, ) 
m lu b D. 

= ''k^k-k-k, aZ-^k-kj-k,, &Y 
2368 367 m 


(2,B,57) 


k^k 2k^kg aW+k^g 3u 


(2 ,B. 58) 
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Equations (2.B.50), (2,B,56) and (2,B.57) are the states of 
interest and are grouped below, 

XI = X2 (2.B.59) 

X2 =-k2k2k^^kgXl"k2]<^^k^:X3^k^k2k2k^^^W (2.B.60) 

X3 = -k2l<3kgkgXl-k3kgk^X3-<k^k2k3kg SW+k^gSu (2.B.61) 

with the output acceleration; 



y = -'k2k3k^^kgXl-k2k^k^X3-k3^k2k3k^ 3W 
Thus the linear system matrices are; 



A 



0 



1 0 



-kjkjk^kj 



-k2k3kgk3 



0 -kgk^k,^ 

0 -Jcjkjk, 



(2.B. 62) 



(2.B. 63) 



B = 



0 

0 

1 



(2.B.64) F=-k3_k2kgk^ 



(2.B. 65) 



G = 



~k]_k2k3k^ 



-k, k„kok_ 
12 3 6 



(2,B.66) 
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Computer simulation of the linear system was accom- 
plished using the Naval Postgraduate School program "lODE' 
operated from the CPCMS terminal of the IBM-360. See Appendix 
C for lODE listing. The program ocean generator for 9W inputs 
a single frequency at a time using the special function fea- 
ture, Wavelength (L,,) will change with encounter frequency 
(co^) as shown in equation (2.B.35), It should be remembered 
that the craft speed is constant at 19,3 knots for the data 
base of this design, 

3 . The Characteristic Equation 

The characteristic equation under no control condi- 
tions of operation, is formed from the A matrix (see Appendix 
A), 



DET(SI-A) 



(2.B.67) 



DET 



s -10 

2267.782 s 11.248 

13002.8 0 64.493s 



(2.B.68) 



Then the characteristic equation is ; 



s^+64,493s^+2267 ,782s+. 5701 (2.B.69) 

which factors into 

(s + 2 , 5414x10''^) Cs+32 . 25+J3 5. 04) Cs+ 3 2 . 2 5- j 3 5 . 04) 

C2.B.70) 
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and is shown in the poler^zero map of Fig. 3. Note that the 
dominant pole is a real pole close to the origin. Also note 
that the complex pair will affect early time response since 
they are relatively far removed from the origin. 

A root locus plot with no control input is shown 
in Fig. 4. The plot covers a range of 

5 < kg < 26 (2.B.71) 

Table II lists the characteristic equation roots for the 
range 



5 < <50 

o 



(2 ,B. 72) 



As can be seen from equation (2.B.63), kg is present 
in both the a^^ and a^g system coefficients. In order to 
perform the root locus study, the characteristic equation 
was rearranged to show the adjustable parameter kg', 

s^+k3k7s4k2k3k^kjSt)<6(k3k^k7k2k3kg+kjk3k^kgk3l<7)tkgS^ 

(2.B.73a) 

with parameter values 



s^+(k.+10 . 5016)s^t2267 ,782s+ ,0928kc 
6 - ■ 0 



(2 ,B.73b) 



From the table and root locus plot, it can be seen that as kg 
becomes large, the complex root pair become real with one 
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root moving toward the origin and the other migrating to 
negative infinity. Additionally, the pole ne^r the origin 
moves very little and remains real throughout the excursion 
of kg . It was observed that this was very similar to the 
root locus plot in the air flow study conducted in Ref, 1. 
However, the removal of the buoyancy term in this study 
caused the pole near the origin to stay real whereas the 
same pole in Ref , 1 becaime immaginary for large values of kg . 
The root locus demonstrated the high sensitivity of the com- 
plex pair to changes in air flow rate, while the pole near 
the origin exhibited a very low sensitivity to changes in 
kg. This observation further demonstrated that the vertical 
oscillatory characteristics of the ship depend upon changes 
in kg. 

As a further check on the effect of air flow on 

C.G. acceleration, an investigation of the effect of changing 

kg on the frequency spectrum of C.G. accleration was conducted 

and the results are shown in Appendix D. From these results 

it seemed appropriate to seek a control law that would adjust 

k_ , the air flow rate. In an effort to gain a better under- 
0 

standing of the effect of various air flow rates on the sys- 
tem, a sensitivity and stability analysis was performed in 
section IVB, The ROOTLO program used in the above study is 
shown in Appendix C. 
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1 




(-32,25 + j35.04) 




2 . 5x10 



-X- 

-4 



X (-32.25-j35.04) 



Figure 3 

POLE/ZERO PLOT FOR THE CI^RACTERISTIC EQUATION- WITH NONI..AL 
AIR FLOW RATE s'^■^64.493s + 2 2 67 . 7 8 2s-<- . 5 7 01 
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X-SCRLE=2. OOE+01 UNITS INCH. 

T-SCRLE = 2. OOE + 01 UNITS INCH-. 

ROOT LOCUS PLOT FOR MODIFIED 
PRESSURE ONLY LINERR MODEL 



Figure 4 
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ONLY MODEL, NO CONTROLS 



I 






I 

i 




ii.i, verification of the model 



A, LINEAR and NONLINEAR SYSTEM VALIDATION 

The validation of the simplified model of a single vari- 
ble pressure lift force was conducted by comparison of its 
frequency response characteristics with those of the model 
presented in Ref. 1, Several techniques were used to 
accomplish this, 

1. Single frequency constant amplitude sinusoids were 
applied to both the nonlinear system of Fig. 1 and 
also the linear system of Fig. 2. The acceleration 
output signal for the nonlinear system is shown in 
Fig. 5. Its sinusoidal form, in steady state indi- 
cates the linearity of the system, which was the case 
throughout the frequency range. 

2. The linear transfer function of center of gravity 
acceleration (Eq. 2.B.15) was subjected to a constant 
wave height input using a spectrum of discrete en- 
counter frequencies. The heave accelerations were 
converted to dB and plotted in a bode diagram shown 
in Fig. 6 

3. Discrete fourier transform technique, developed in 
an independent study by S, Carpenter, [Ref, 9], was 
also used to test the response of the linear system 
(Fig. 9) to an irregular sea input shown in Fig. 7 
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consisting of twelve components of frequency and 
amplitude from the Pierson-^Moskowitz energy density 
spectrum given in Fig. 8 [Ref. 8, 9], 

The results of procedures 2 ^nd 3 above were compared 
with the linear system response of Ref. 1, The magnitudes 
for the center of gravity a,cceleration showed approximately 
20% difference between the pressure lift only, no control, 
response curve (Fig. 6) and the linear system of Ref. 1 
(Fig. 10). The general shape of the curves are identical. 
Therefore, the pressure lift only model was assumed to be 
valid and useful for linear regulator control design. The 
observed peaks and nulls of Figs. 6, 10 are a result of the 
sine function numerator coefficient, k^, as described in 
Section III. The relationship between plenum length 
and wavelength (L ) expressed in k, , determine where, in 
the encounter frequency spectrum, the nulls and peaks occur. 
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SYSTEM STRBIUTT TESTS 
oBODE PLOT, MflGIDB) VS, LOG N 
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Figure 6 
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Figure 7 
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Figure 8 
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IV, THE CONTROL LAW- 



A. LINEAR REGULATOR DESIGN APPROACH 

The linear pressure lift only model is described in 
state equation form; 



— — 








— 








— 
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0 








0 
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^21 


0 


^23 






+ 


0 






^31 


0 


^33 








1 


— 





















with the cost function of the form 



J = 




(x'^QX + 9U'^§3U) 



(4. A. 2) 



where the final value error is not included. The system may 
be designed using the linear regulator technique with the 
given quadratic cost function. An optimal control was 
found using both the transient and steady state solutions 
of the Riccati equation. The control law 



" .^1 T 

8U =?R B KX 

-V -- — 



C4.A, 3) 



where the K matrix was found from the solution of the Riccati 
equation 
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4 



^ ~ - Q + K B 

" ^ n r* 



(4, A. 4) 



with the appropriate choice of the Q and R coefficients, the 

n rt 



feedback gains for the system were calculated. Unfortunately, 
the first choice for Q and R may not yield control gains which 



are physically realizable or practical. For this reason, it 
is appropriate to examine the effects of the control gains on 
the pole location in the s plane for the controlled system. 
Such a study was done and is presented in part B of this 
section. The choices for Q and R were derived by trial and 
error . However y 1:he design crilierion was "to reduce this mag- 
nitude of the frequency response curves through the study of 
various Bode plots. To this end, many Q and R values were 
tried . 

The Riccati equation computations were done using a com- 
puter program documented in Ref, 10 and modified by B. 
DesJardins in Ref. 11. Appendix C illustrates the input/ 
output format of this program. Based upon the above design 
criterion, the final weighting matrix chosen was j 



10 



5 



0 



0 



Q 



0 



10 



3 



0 



(4. A, 5) 



0 



0 



. 5 



and 



R = . 03 
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The optimal control ga,ins obtained using the above Q and R 
matrixes were; 

= -'1,689x10^ 

g^2 = 1.19x10^ (4. A. 6) 

g ^3 = -1.83x10^ 

These gains were multiplied by their respective states and 
applied to the state equations (2.B.59, 61, 62). Under state 
feedback control, the resultant output acceleration frequency 
magnitude was reduced by 13% from the peak no control magni- 
tude. Table III lists the peak lobe values of the C.G, 
Acceleration/wave height magnitudes versus frequency for the 
nonlinear and linear systems. Each system was exercised 
under control and no control conditions. The steady state 
sinusoidal output for the linear model, under control and 
no control conditions is shown in Table IV for (o^ = 45 rad/sec, 
where 



^u “ ^11"^ ^12^2 ^13^3 

and 

y = *3Z 

The control signal for the nonlinear system was generated by 
using the linear control low coefficient found using linear 
regulator theory. The reduction in acceleration was not 
uniform throughout the frequency range. At law frequencies, 

(o^<17 rad/sec, control gains caused an increase in C.G. acceleration. 



(4. A. 7) 
(4. A. 8) 
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As 0 )^ APP^i^oached 100 ra.d/seCj C,G, acceleration approached but 
never exceeded the no control condition. Effective control 
was realized in the frequency range 17<W^<100 radians/ sec, 

B, SENSITIVITY AND STABILITY 

Since a set of optimal gains were found for the system, 
it remained to test the system under feedback control and 
observe the sensitivity and stability for acceptability. 

The locus of roots of the characteristic equation for 
the pressure lift only system without state feedback control 
was shown in Section III, The objective of state feedback 
control is to adjust the position of the complex pair of 
roots in such a way as to reduce center of gravity accelera- 
tions. For the open control loop, the transfer function 
was shown (Eq. 2,B,22) to be 



3W 



2 

-as 



s +bs+c 



(4.B.1) 



Addition of state feedback will change the value of b and c. 
Note that if the k^ sine function null effects are removed 
from the numerator , then a=c . It follows that the magnitude 
of the Bode plot will reach unity value at high frequency 
regardless of changes in b and c. In other words, the sys- 
tem has high pass frequency characteristics. Thus, the 
results of increasing the value of c will be a reduction of 
the magnitude plot in the lower range of frequencies of 
interest. Also, note that changes in the b coefficient will 
affect the resonant peaking at the corner frequency, 
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TABLE III 

PEAK LOBE VALUES OF FREQUENCY RESPONSE, PRESSURE LIST ONLY 
LINEAF S NONLINEAR MODELS UNDER CONTROL AND NO CONTROL 

CONDITIONS 



Wg NONLINEAR LINEAR 





NO CONTROL 


CONTROL 


NO CONTROL 


CONTROL 


6 


18,6 


26,5 


18.5 


27 , 8 


7 


18,5 


25,2 


18 . 6 


26.6 


8 


15.7 


20,3 


15.7 


21.4 


9 


9.7 


12.7 


9.8 


12 . 7 


14 


42.2 


42,2 


42.3 


45 . 7 


15 


48.2 


48.2 


48 . 4 


50 . 9 


16 


49.7 


49.7 


49,9 


51.1 


17 


45.9 


45 , 9 


46.1 


46 . 2 


24 


69.8 


62.9 


70.5 


63 . 4 


25 


78.1 


69.7 


79.0 


70.3 


26 


79.1 


70.0 


79 . 9 


70 . 6 


27 


72.1 


63.4 


72 . 8 


63 . 9 


28 










34 


90.8 


77.7 


91.8 


78 . 5 


35 


100 . 8 


86.3 


101.8 


87.3 


36 


102.4 


86 . 8 


102.3 


87.9 


37 


91.7 


78.6 


92 . 6 


79.6 


38 










44 


99.5 


86 , 8 


100.1 


87.4 


45 


110.0 


96.4 


113.0 


98 . 6 


46 


110.0 


96 , 8 


111.5 


98,0 


47 


99.7 


88 , 0 


101.0 


88 . 9 


48 










54 


97.1 


88 . 5 


97 . 5 


88.5 


55 


107,7 


98,5 


108 , 0 


98.7 


56 


108.1 


99 , 2 


108 . 8 


100 , 0 


57 


98.2 


90 . 5 


99.0 


91.2 
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0) 

e 

64 

65 

66 

67 

68 

69 

74 

75 

76 

77 

78 

79 

84 

85 

86 

87 

94 

95 

96 

97 



TABLE III 
(continued) 



nonlinear linear 

NO control control no control control 



88.7 


83,8 


89,3 


84 . 3 


99.3 


94,1 


99.5 


94 . 8 


100.4 


95,4 


101.1 


96 , 5 


91,9 


87.6 


91.9 


88 , 0 



79 . 0 


76,5 


79,4 


77,0 


89.5 


86 , 9 


89 . 9 


87,3 


91,3 


88,8 


91.3 


88.6 


84.4 


82 , 2 


84.8 


82.6 



69,9 


68.8 


70.1 


69 , 0 


80 , 2 


78.9 


80.1 


79,3 


82.6 


81.4 


83.0 


81,7 


77,4 


76.1 


77.0 


76,5 


61.9 


61.4 


62.0 


61.5 


71.9 


71,5 


72.1 


71.6 


75 , 0 


74.5 


75.2 


74.7 


70.9 


70.5 


70.9 


70.6 
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TABLE IV 

SINUSOIDAL LINEAR MODEL OUTPUT WITH/WITHOUT FEEDBACK CONTROL 

AT ci)^=:45,0 Rad/Sec. 



TIME 

T 


CONTROL 

Y 


CONTROL EFFORT 
3 

u 


NO CONTROL 
Y 


0.780 


54.38 


4,28 


61.02 


0.785 


34.33 


2 . 64 


38,03 


0.790 


12 , 50 


0.86 


13 . 08 


0.795 


-9 . 96 


-0,95 


-12,54 


0.800 


-31.91 


-2,7 2 


-37.52 


0.805 


-52.23 


-4,35 


-60.56 


0 . 810 


-69,85 


-'5 .76 


-80.49 


0.815 


•^83,86 


-6.87 


-96.26 


0,820 


-93.56 


-7.62 


-107 . 08 


0.825 


-98.44 


-7,99 


-112 ,38 


0.830 


-98,23 


-7.94 


-111.88 


0.835 


-92.97 


-7,48 


-105.62 


0.840 


-82.90 


-6.64 


-93 .92 


0.845 


-68 . 60 


-5,45 


-77.37 


0.850 


-50.74 


-3.97 


-56.85 


0.855 


-30 , 27 


-2.32 


-33.39 


0.860 


-8 . 34 


-0.53 


-8 . 21 


0.865 


14.22 


1.29 


17.40 


0 . 870 


35.95 


3.04 


42.10 


0,875 


55.82 


4.63 


64.64 


0,880 


72.82 


5.99 


83.85 


0.885 


86 , 06 


7.03 


98 .73 


0.890 


94.87 


7 .72 


108 . 53 


0.895 


98.79 


8.00 


112.73 


0,900 


97.61 


7,88 


111.11 


0.905 


91,42 


7.34 


103.79 


0,910 


80 , 50 


6.43 


91.11 


0.915 


65.44 


5.19 


73.74 


0.920 


47.01 


3.68 


52 . 56 


0,925 


26.49 


1.98 


28 . 68 


0,930 


3,94 


0,17 


3 . 31 
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1 • The Characteristic Equation For Feedback Control 
A pplie d 

The signa,! flow graph of Fig. 2 was modified to show 
feedback gains a,nd is shown in Fig, 11. Using Masons 
gain rule, the characteristic equation was derived from the 
S.F.G. (Eq. 4,B.10). The resulting equation was compared 
to characteristic equation derived from the modified system 
matrix with feedback of Appendix D. It was noted that they 
differ by the terms 



^ 21^3 3 '"^2 3^31 '■ 



These terms very nearly cancel and are insignificant (a2i 
a 33 ~a 23 a 3 i= . 5701) when compared with the other larger terms 
in the characteristic equation. Therefore, their omission 
from the signal flow graph of Fig. 11 is valid. The deri- 
vation follows: 

. _ -^3 



”®3 _ ~^3 



(4.B.2) 



(k3)(-kg)(-k.7) 



"^3^6^7 



-a 
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(4.B. 3) 



(-k2)(k3)(-k^^)(-kg) 



”^2^3^4^8 



^3. 



(4.B.4) 



(k.^)(k3)(-k^^)Cg2) 






_ ^23 



f g2 (4,B.5) 



(k^) Ckg) (-k^^) 



^3^4^7 



©1 



'23 



f gl C4.B,6) 
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Figure 11 

Linear Pressure Lift Only 
Signal Flow Graph With State Feedback Control 



A = I-CL^+Lj+Lj L3 



(4.B.7) 




a 



s 



23 



3 




s 



(4 .B. 8) 





(4.B.9) 




(4 .B. 10) 



2 , Root Locus 

The following root locus plots are a study of the 
effects of changing the three different control feedback gains. 
Many weighting combinations were tried and it was discovered 
that the system was most sensitive to changes in air flow rates. 
This observation was in agreement with the analysis of Section 
II that the system was most sensitive to kg . Each of the 
following root locus plots and associated table of values was 
chosen as representative of the alternatives tried. The 
characteristic equation used is shown as Eq, (4.B.10). The 
Naval Postgraduate School program for the IBM 360/67 computer 
was used to generate the tables and plots. The program 
listing is shown in Appendix C, 
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In the first study ^ the constant gains used were 



~ 425,0 j ~ 

and veried over three decades starting from g „ = .1. 

The resultant plot and table of values are shown in Fig. 12 
and Table IV respectively. For ease of reference, the syS' 
tern poles were given letter designations for the characteris- 
tic equation 

(S+A) (S+B) (S+C) = 0 (4.B.11) 

where 

A = (+32.25-j35.04) 

B = (+32.25+j35.04) 

C = (2. SxlO-'^ + jO) 

As g ^2 nears -20, pole a rapidly enters the real axis and 
migrates towards negative infinity and pole B moves upward 
across the real axis, Pole C moves in a negative direction 
until g ^2 = *'26 and then moves below the. real axis tracking the 
mirror image path of pole B. It was seen from this plot 
that as g ^2 nears -20 that the complex pair undergoes the 
greatest movement, In this study there was no point at 
which instability occurred, When the gains were tested on 
the syste.m, there was no notic^ble decrease in vertical 
acceleration. 

The next phase of this investigation uses constant 

gains of 
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X-SCflLE=2. OOE+01 UNITS INCH. 
T-SCflLE=2. OOE+01 UNITS INCH. 



ROOT LOCUS PLOT FOR MODIFIED 
PRESSURE ONLY LINERR MODEL 



Figure 12 
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TABLE IV 



Roots For 


Modified Pressure 


Lift Only Model , 


gll-425.0, 




®13 


Ya,ries 




A 


B 


C 


-0,099 


-31.49+j46 .96 


-31.49-.j46.96 


-1 . 4 9 + j 0 


-0,127 


-31,52+j46.88 


-31.52-j46.88 


-1. 586 + j 0 


-0.162 


-31.52+j46 ,86 


-31.52-j46.86 


-1. 611+jO 


-0.207 


-31.53+j46.84 


.-31,53-j46.84 


-1. 643 + j 0 


-0 .263 


-31,53+j46,81 


-31. 53-j46 . 81 


-1 . 684 + j 0 


-0,335 


-31,54+j46,76 


-31, 54..j46,76 


-1. 736 + j 0 


-0.427 


-31.55+j46,71 


-31.55r,j46,71 


-1. 802+j0 


-0 . 544 


-31. 57+j46 ,64 


-31, 57-j46 . 64 


-1. 887+jO 


-0 , 693 


-31. 59+j46 , 56 


-31. 59-j46 . 56 


-1.995+jO 


-0 , 883 


-31. 61+j46 .45 


-31.61-j46.45 


-2 .134+jO 


-1.125 


-31.64+j46 .30 


-31,64-j46.30 


-2.312+jO 


-1.433 


-31,68+j46.12 


-31.68-j46.12 


-2.540+j0 


-1.825 


-31.72_j45.88 


-31,72-j45.88 


-2,835+jO 


-2,326 


-31.77+j45. 58 


-31,77-j45.58 


-3,216+jO 


-2.962 


-31.84+j45.18 


-31. 84-j45 .18 


-3.711+jO 


-3.773 


-31. 91_j44 .67 


-31,91-j44.67 


-4 . 359 + jO 


-4.807 


-31,99+j43.99 


-31 . 99-j43 .99 


-5.214+jO 


-6.124 


-32.06+j43,08 


-32 . 06-j43 . 08 


-6,359+jO 


-7.801 


-32.09+j41.87 


-32 . 09-j41. 87 


-7 . 923 + jO 


-9.937 


-32.03+j40.20 


-32.03-j40.20 


-10.134+j0 


-12.658 


-31.70+j37.87 


-31.70-j37,87 


-13 .450+j0 


-16.125 


-30 . 60 + j 34 . 50 


-30.60-j34.50 


-19.032+j0 


-20 . 540 


-26 . 95+j30 . 22 


-26.95-j30.22 


-30,645+j0 


-26,165 


-48.97+jO 


-20.53+j29.31 


-20.534-j29. 31 


-33,331 


-64. 98+jO 


-16.02+j30,86 


-16.023-j30.86 


-42.458 


-80.26+j0 


-12 . 84 + j 32 .65 


-12 . 838-j 32 ,65 


-54.085 


-96 . 55+jO 


-10 , 36 + j 34 . 38 


-10.367-j34,38 


-68 . 90 


-115,02+j0 


-8 ,36+36 . 02 


-8 . 363-j36 . 02 
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34,98 



gll = -195,32, 



§12 



while varied over two decades from -.1. These gains 

were calculated from the Riccati with weighting matrix*, 



Q 



lo"^ 0 0 

0 12 ^ 0 

0 0 .5 



and 



(4.B.12) 



R = ,03 



(4.B.13) 



The resulting plot and table of pole values are shown in 
Fig. 13 and Table V respectively. As varies, the poles 

follow similar paths as in the previous study. The greatest 
degree of pole movement, or sensitivity, occurs at g^^""!®**^* 
This set of gains produced instability when 

In the final study, the constant gains were 



g^^ = -1688.095, g^2 = 119.16, g^3 = -18.306 

Again, g^^^ was made to vary over three decades starting with 
g ^^3 = -.1. The resultant plot and table of values are shown 
in Fig. 14 and Table VI respectively. These gains were cal- 
culated using the Q and R values shown in (4,A.5). Again, 

'1 

the poles follow similar paths but pass closer to the origin 
than in the previous two studies, causing a greater net 
decrease in the systems oscillatory tendencies, As in the 




X-SCRLE=2. OOE+01 UNITS INCH. 
T-SCflLE=2. OOE+01 UNITS INCH. 



ROOT LOCUS PLOT FOR MODIFIED 
PRESSURE ONLY LINERR MODEL 



Figure 13 " 
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TABLE V 



Roots For 


Modified Pressure 


Lift Only Model , 


g^^ = -195. 32 j 




g^2=34 ,98 


, g^ 3 =varies 




A 


B 


C 


-0.099 


-32,65+j40,60 


-32.56-j40,60 


0, 81+ jo 


-0.127 


-32 , 66 + j40 , 50 


-32,66-j40,50 


0.71+j0 


-0,162 


-32 . 67+j40 ,48 


-32. 67-j40,48 


0. 68 + jO 


-0,207 


-32,67+j40,45 


-32,69-j40,45 


0.64+j0 


-0,263 


-32. 67 + j40 .41 


-32 . 67-j40 ,41 


0. 60 + j0 


-0.335 


-32,68+j40.35 


-32 . 68-j40 .35 


0. 54 + jO 


-0.427 


-32.69+j40.29 


-32.69-j40.29 


0.46+j0 


-0.544 


-32,70+j40.20 


-32.70-j40.20 


0.37+j0 


-0.693 


-32,71+j40.09 


-32.71-j40.10 


0 , 25 + jO 


-0.883 


-32,72+j39.95 


-32.72-j39.95 


-.09+j0 


-1.125 


-32.74+j39,78 


-32,74-j39,78 


-Q.ll+jO 


-1.433 


-32.76+j39.54 


-32 ,76-j39. 54 


-0. 37+jO 


-1.826 


-32 .78 + j39. 25 


-32.78-j39.24 


-0. 71+jO 


-2 .326 


-32. 81+j38 ,86 


-32 ,81-j38 . 86 


-1.14+jO 


-2 , 962 


-32,84+j38.35 


-32.84-j38.35 


-1. 71+jO 


-3 . 774 


-32.86+j37.68 


-32,86-j37.68 


-2.46+jO 


-4 .807 


-32 . 86 + j36 .79 


-32,86-j36,79 


-3 ,47 + jO 


-6.124 


-32 , 81+j35 . 59 


-32,81-j35,59 


-4 . 8 5 + j 0 


-7.801 


-32,66+j33,93 


-32.66-j33.93 


-6.80+j0 


-9 . 937 


-32 . 23+j31, 54 


-32.23-j31.54 


-9.74+jO 


-12.66 


-31.01+j27.93 


-31.01-j27.93 


-11.48+jO 


-16.12 


-26 . 18+j22. 73 


-26.18-j22.73 


-27 . 87+jO 


-20 . 54 


-47,81+jO 


-l8.37-j23.83 


-18 . 37-j 23 . 83 


-26.16 


-61.49+jO 


-14,27-j26.51 


-14.27-j26.51 


-33,33 


-74 ,16+jO 


-ll,43-j28. 89 


-11. 43-j 28 ,88 


-42,46 


-87,51+jO 


-92,12-j31,05 


-9,21-j31, 05 


-54.09 


-102,51+j0 


-73,88-j33.04 


-7,39-j33,04 


-68.90 


-120.00+j0 


-58.68-j34,90 


-5.87-j34,90 
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X-SCflLE=2. 00E»01 UNITS INCH. 
T-SCflLE=2.00E+01 UNITS INCH. 



ROOT LOCUS PLOT FOR MODIFIED 
PRESSURE ONLY LINEAR MODEL 



Figure 14 
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TABLE VI 

Roots For Modified Pressure Lift Only Model, gT„=119, 

t^^Vdries) 



Hi, 


A 


B 


C 


-0.099 


-34. 66+j52 ,35 


-34,66-j52,35 


4,817+jO.O 


-0.127 


-34.68+j52.29 


-34.68-j52.29 


4.751+jO.O 


-0.162 


-34,69 + j52 ,27 


-34.69-j52,27 


4,733+jO. 0 


-0 .207 


-34.70+j52.25 


-34,70-j52,25 


4,709+j0.0 


-0.263 


-34,72+j52,22 


-34,72-j52.22 


4.680+j0.0 


-0 .335 


-34.73+j52.19 


-34’,73-j52.19 


4.642+jO.O 


-0.427 


-34.75+j52.15 


-34.75-j52,15 


4.594+jO.O 


-0 . 544 


-34.78+j52,09 


-34.78-j52,01 


4.533+jO.O 


-0 . 693 


-34.81+j52.02 


-34.81-j52.02 


4,454+jO.O 


-0,883 


-34.85+j51.93 


-34. 85-j51,93 


4 . 355 + jO . 0 


-1.125 


-34. 91+j51, 81 


-34.91-j51. 81 


4,227+jO.O 


-1.433 


-34,98+j51.66 


-34.98-j51.66 


4,064+j0.0 


-1.827 


-35,06+j51.46 


-35.06-j51.46 


3.854+jO.O 


-2 .326 


-35.17+j51,21 


-35,17-j51.21 


3 . 585 + j 0 . 0 


-2 . 962 


-35.31+j50.89 


-35.31-j50.89 


3,240+j0.0 


-3 . 773 


-35.48+j50.46 


-35.49-j50.46 


2,794+jO.O 


-4.807 


-35.70+j49.91 


-35.70-j49,91 


2.216+jO.O 


-6.124 


-35.97+j49.19 


-35.97-j49,19 


1.463+jO.O 


-7.801 


-36.29+j48.22 


-36 , 29-j48 .22 


0.472+j0.0 


-9 . 937 


-36.67+j46.91 


-36.67-j46.91 


-0.849+j0.0 


-12 . 66 


-37.10+j45.11 


-37.10-j45.ll 


-2.648+jO.O 


-16.12 


-37 . 52+j42 .55 


-37,52-j42.55 


-5.191+jO.O 


-20.54 


-37.74+j38,72 


-37.74-j38,72 


-9. 056 + j0 . 0 


-26.16 


-36.96+j32.14 


-36,96-j32,41 


-16.11+jO.O 


-33,33 


-42,45+j0,0 


-27,29+j23 .37 


-27.29-j23,37 


-42.46 


-69 , 58 + jO, 0 


-18 ,18+j27 .34 


-18,18-j27.34 


-54,09 


-8 9 , 5 6 + j 0 , 0 


-13.86+j30,54 


-13.86-j30. 54 


-68 . 90 


-109.98+j0.0 


-10.88+j33.10 


-10.88-j33.10 
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previous study, some values of caused instability. 



Particularly, for pole C was in the right half 

plane, The greatest sensitivity was experienced for values 
of g ^2 riear g^^ ~ -SO. Note that the value of g^^ chosen 
via Riccati will cause stable operation and is also removed 
from the value of greatest sensitivity. 

3 . Routh Criterit>n 



the range of feedback gains for stable operation obtained 
from the root locus study. The Routh array is set up below. 
The C.E,, given in Eq. 4.B.17, is used with the understanding 
that all coefficients are absolute values. 



The Routh criterion was applied in order to confirm 



1 



^21 ^23 ®2 



^33 ®3 



^"^23 ®1 * ^21 ® 3 ^ 




0 



~^ 23®1 * ^ 21^3 



0 



The conditions for stability are: 



1) gg >. ^a33 



(4.B.14) 



since a 33 - 64,493, then 
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(4.B.15) 



^ 21^3 ^ ^ 23^1 

a . 



Si < 



"21 



'^ 23 ' ®3 



^ 2267,782 
^1 11/248' ' §3 



< 201,616 g3 

g^ < (201. 616) (-64.493) 

g^ < -13002,821 

. X . ... ,^ 23 Sl ''^ 21 ^ 3 , 

3) 32^ 323 g2 > ”( = rt:^ ) 



^ 23 '^ S 3 



So > 



. ^ 23 ^ l '^ 21 ^ 3 v 

I 



- a 



21 



'23 



^ ^ 23 ^ 1 ~^ 21 ^ 3 ^ ^ ^ 21 ^ 33 ~^^ 21 ^ 3 ^ 



g 2 > 



^ 23'^®3 



^ 33-^23 



23 



§2 " 



'^ 23 ^ 1 ~^^ 21 ^ 3 ~^ 21 ^ 33 ~^ 21^3 

^ 33~*~^3 

^23 



So > 



"^33^l"^21^33 

^ 33'*’^3 



"23 



^^ 33 ^ 1 ~^^ 21 ^ 33 ^ 

^ 23 ^^ 33^®3 



(4.B.16) 
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If gg > -64 the g 2 > -201,72, 

The bounda,ry limits of g^, and therefore g^ , are 
very sensitive to g^. 

For the specific control chosen in Section IVA, the 
Routh array was used to determine the stability limit for 
the gains. 

From Eq(4,B.14) 

®13 ^ 



therefore 18.3 > -64,493 

From Eq( 4 . B , 15 ) 

. ^21 ^ 

§11 §13 

therefore 1688 < 3689,581 



From Eq(4.B.16) 



^33§ll'*‘^21^33 

^ 23 ^^ 33''’§13 



therefore 119 > -273,95 

It was seen from this study that the gains were well 
within the limits of stability. Lower gains were realitively 
more stable. 
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V. CONCLUSION 



It has been shown that the linear regulator design can 
be applied to the problem of vertical plane motion for the 
captured air bubble surface effect ship. The linear pressure 
lift only model, developed from the simplified equations of 
motion, was shown to be a valid model which could be used for 
control analysis and design. An optimal control was found, 
based upon a quadratic cost function, using empirically de- 
rived weighting factors , The pressure lift only linear 
system showed a reduction in center of gravity acceleration 

for 17 <o) <100 radians/sec, A C.G. acceleration reduction 
e 

was observed in both the linear and nonlinear systems . The 
shape of the frequency response magnitude curves were identi- 
cal. A 13% maximum reduction was observed at oo =45 rad/sec. 

e 

where peak acceleration occurs . 

The control gains derived, as shown in Section IV, are 
within the airflow rates currently available for the XR-3. 

The Riccati equation computer program on Ref. 10 proved to 
be a most useful tool, which facilitated experimentation of 
a large number of weighting coefficients , making an optimal 
control design more feasible, 

Finally, the analysis shows that the system is stable 
over the range for the control gains selected. It was seen 
from the root locus study that the limits of stability narrowed, 
that the sensitivity to air flow rates increased and C.G. 
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h 



acceleration decreased with higher control gains, It was 
seen that for these higher gains the C,G, acceleration was 
not reduced and the control, 3 , increased for u <17 rad/sec. 
This was due to the fact that the buoyancy effects are pre- 
dominant in that frequency range and the buoyancy terms are 
at a constant steady state value in this model. 

Though an optimal control was found based on a particu- 
lar quadratic cost function, many other possibilities exist. 
Different defined weighting factors can be experimented with, 
and some may result in improved C,G, reduction. The limiting 
factor in such a design is equipment limitation. The choice 
of weighting factors, and ultimately, the control, are even- 
tually determined by that fact. 
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VI. RECOMMENDATIONS FOR FURTHER STUDY 



The following list of reconutendations is presented for 
consideration and further study. 

A. A more realistic system presentation could be accom- 
plished if pitch coupling was combined into the 
linear regulator study. This would facilitate the 
speed analysis study for the avoidance of wet decks 
and undesirable pitch accelerations. 

B. The passive control technique discussed by Boggio 

in Ref. 6 should be examined further. The introduc- 
tion of a flexible membrane into the plenum chamber 
may reduce vertical accelerations. This technique 
has been examined in model towing tank studies, 

C. Since the mass of plenum air X3 cannot be directly 
measured, an observer might be used in conjunction 
with state X3 or perhaps remodel using low order 
model analysis. Perhaps consideration should be 
given to a design using only XI and X2 as states to 
be used for feedback, 

D. The linear regulator technique design approach should 
be compared in greater detail to the filter design 
technique proposed by Rohr Marine Inc, in Ref. 3, A 
fixed air flow rate, based upon a filtered observed 
sea state, could be a more realizable design in terms 



67 



of the physical equipment limitations* The study 
of fixed air flow, of Appendix E, shows that the 
filter design approach may offer greater control 
over the entire encounter frequency range, (see 
Fig. E-1) 
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APPENDIX A 



The simplified non linear system of Ref. 1 is shown in 
Fig, A1 . This system, with the buoyancy term removed, served 
^s the basis for study in this thesis. The linearization 
and resultant linea,r system matrix is shown for the simpli- 
fied non linear system, The pressure lift only model lineari- 
zation is also shown. 

LINEARIZATION OF THE SIMPLIFIED MODEL 

The equations of motion are linearized using Taylor 
series expansion around the steady state operating point. 

Heave acceleration: 
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In a similar manner 




figure A1 



70 



8 q. = '(n-^) 3 'P, 

in kq b 



3 q = . ,p 

^ 2F,C0) t' 

D 



~ vr(o) ^^b^ 

m b 



'''b = 



3P, = 3P. 
b b 



31 h = 9 , 

d z 
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yp ( 0 ). 
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r 51 b T 

'■T^coy " Tjyy^ 

in b 



3V = -[n +^i yP (0), 

"" 2 F^( 0 ) ^ 

b 

m b 

V~TF)' 

m b 

The state vector is 

X = [3Z 3Z 3V ]'^ 

— m 

The output variable is 
y = 3P^ 

^ ~ ^ 11^1 '*’ ^ 13^3 

where 

yP, (0)A, 

r _ b b 

^11 VJb)- ■ 

b 

,P-( 0 )A ^2 

^13 r~TW~^ 
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Writing the above equations in state equation form; 



X 



1 



= X 



2 



X 



2 



2A P 1 

r r s w' 

ir- 



+ 



YP, (0)A, ^ 

( ^ ) 1 X 



yP, (0)A, 
m 



• ^ P 

X. = [-YP>.(0)Cnk ^ 



^a 2P7(0) 

D 






P ' ~l y 

V^roT-* 



* t-YPt,(0)CnY^( 



V 1/ 2 -| ^ X 

2P.(0) ! 



or 



^2 ■ ^ 21^1 ^ 23^3 

^3 ■ ^ 31^1 ^ 33^3 



The characteristic equation was found from the system mi^trix , 



DETCsI-A) = 0 
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to be 



+ 62s + 2183s + 789 = 0 



which factors into j 



(s+.37)(s+30.82 + j34,8)(s + 30.82T^j34,8) 



TAYLOR SERIES EXPANSION OF PRESSURE LIFT ONLY VERTICAL 
EQUATIONS OF MOTION 



Defining 





F 



pm 



F 

P 

M 



then incremental acceleration due to pressure is . 



Z (0) + dZ 



= W - (F (0) 
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.. 9F 
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Defining the air mass volume: 
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Then the incremental volumetric change in air mass is 



8V 

m 



(nk )9P? 
q b 



^n^l ^ . P'a )l/2 

^a 2PT(0) 
b 




and the incremental force due to pressure isj 



8F 



P 



b b 



Incremental air flow into the plenum 




- ( nk ) SFr 
q b 



Incremental air flow out of the plenum 
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Incremental gauge pressure (plenum) 
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Incremental plenum volume 



9V. 
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Incremental Pressure (plunum) 




3P 



b 



Incremental Draft 
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Acceleration (3Z) 
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Volumetric air mass ( 3 V ) 
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Plenum Pressure 



9V 



3Pb = yP^CO) 



m 



A series of system coefficients , k , were defined for ease 



n 



of writing and are used throughout the thesis. They are. 
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APPENDIX B 



A brief step weight removal study was performed and 
compared to Ref. 1 as additional check on system validity. 
The lODE remote terminal CP-^CMS program on the IBM/ 360 was 
used in this study. A sample listing is shown in Appendix 

C. 
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Figure Bl 
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APPENDIX C 



The following are the listings for the various computer 
programs used in the preparation of this thesis. 



\ 
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OPTIMAL CCM*^OL/K.tLMAM rlLTEP POCGPAM 
PR08LSM lOrINTIFICATICM - HSAVS CINTROL 






THE A MATRIX 



0.0 

-2. 267782100 C3 
-l.3002bC70C C4 



TH: e MATRIX 



l.OOOJOOuOD 00 0.> 

g.O -1.12A300000 01 

C.C •6.4^S300C00 01 



0«0 C «0 

TH^ C MATRIX 

•2.26778200C 03 0.0 



l.OOCQOOQCO 00 
-1.124800000 01 



4 « « H :^c « 4> 4e 4e* 



CONTROL OPTION *** 



THC R MATRIX 

3.C00Q0000C-02 

THE C MATRIX 

l.QOOOOOOOC 05 0.0 

O.C' l.OJOOJCOOC ^>3 

0.0 C.O 

STEADY STATE SOLUTION 

GAINS 

1.68809469C 03 -1.191156:60 02 



O.C 

5.00C00000D-01 



1.830578260 01 






♦♦♦ CONTROL OPTION 



THP R MATRIX 

i.GOOOOOOOC-02 

THt' 0 MATRIX 

l.JOOO^OOOC 05 0.0 

O.C l.COCCCOOOO 03 

0.0 O.Q 

INITIAL CCNOITICNS 



0.0 

0.0 

O.C 



0.0 

0.0 

0.0 



1.0 

0.0 

5.CCC00000D-O1 



G.O 

o.O 

C.O 



3|S^ ^ ^ * ♦ V * # ^ ^ ^ ^ ♦ 



TIME 

GAINS 



0.0 



2.CCC0Q1CCC Cl 

0.0 



TIME = . 
GAINS 

1 .638J95J6C 

TIMc 

GAINS 

1.688J9469C 

TIMC =♦ 
GAINS 

1.638C9469C 

TIME *• 
GAINS 

1.688C9469D 



1.9CC00002C 01 
03 -1.191155710 

1.6CC0COC40 01 
03 -1.19115636C 

1.7CCGC0C6C Cl 
03 -1.191156C6C 

1.6CCCCCC£D Cl 
C3 -1.191156060 



T2 

32 

02 

02 



0.0 

1.830577720 01 
1.830578260 01 
1.830573260 01 
1.83357826C 01 
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The INPUT CATA IS 

GRAPH TITL'= A.‘iO SCALIHG DATA 

ROOT LfCLS PLCT FTP wCCIFI'C 
PRESSUk- only LiN?:AR VCCRL 



XSCALP 


YSCAL‘= 


!XUf 


lYCIGh 


lU ICE 


IhIGh 


IGR to 


20.0 


?0.0 


1 


6 


7 


6 


1 



ORDFR CF TH- Ch/- AC TER I < T IC r^CLAHCN 

3 



constant P^RT The CCEFFICIF-MS IN 03SC^NDTNG CPD5.R 
O.IOCGCO 01 C.T 0.2Z61EC CA C.57010C CO 



VARU3LF FART OF Th- COEFFICIENTS 
0*0 O.l'JACTD 32 C,0 

INITIAL VALLF Cl* ThS V^FIAeiE 
C,500CGE 01 



IN CESC'^'MCrNG CFOIER 

c.c 



NUM8CP OF C=CAjrS TO 9^ SPANMSC 

1 

THE SYSTEM PuL*' aPe 

I N S L(U'i<TST) VsAPN INC (U ITh FIX') 

REAL PART 

3.12 570F-C3 J.125TOS-C3 -0 . 2 5 139 

IrtAGIMARV PART 

O.A762l<^ C2 -C.A7621^= C2 C.C 
THE SYSTEM ZFROS i: 



IQ.MR 



REAL PART 

(T.r> 0.0 

IMAGINARY PART 



2 (I HR = 6o) 



0.0 



C.C 



ALL POINTS HAV- T^F CCC F C I N AT E S . 

THIS PLOT V^ILL NOT 9F ClTFLT. 






RR 



SPECIFICATIONS 

VARIABLES 6 INITIAL CONCITIONS: 

XL = 0,0 
X2 = C.C 
X3 * 0.0 
T ^ C.C 

CONSTANTS: 

A21 = 2267.7fl2CC0 
A23 = 11.2A8COCOO 
Ail = 2llt.245CCC 
A33 = 10.A97C00C0 
B2l = 2267.782000 
621 = 2116.265C00 
832 = l.CCCOCOCOO 
G = 32.17AC0C00 
A = C.ICCCCCCCCOD 00 
WE = A5.5CCCCCCC 
K6 6.1AA22C0C0 
Kll = -i6L6.C9A700 
K12 * 119.I160C00 
K13 = -16.3C6CC000 

SPECIAL FU^CTICNS: 

LW = (2.*3.14159265*GJ/WE 
F * (SIM£2.£319/LWJ)/<62.8319/LW) 

SEA = A*SINIWE*T) 

U2 = KU’»X1+K12*X2+K13*X3 
y = -A21*Xl-A22*X2-B2l»F«SEA 

DERI VATI VES: 

CtXl /C<T ) = = 

X2 

D(X2 /OIT ) = =* 

-A2i*Xl-A23*X3-821*F*SEA 
D(X3 /C(T ) = = 

-A3l^K6^X l-A32*K6*X3-03l*K6*F+SEA+U2 
OLTP UTS: 

TITLE: CAB SES INCREASEC AIR FLOW RATE 
tabulate: 7 Y 

AT interval 0.5COOCOOOOOD-02 
NG FLCTTING 

ENC CALCULATION WHEN T .GE . 1.50000 
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—.—specification—^-^- 

VARIABLES S INITIAL CONDITIONS; 
XI = 0.500000 
X2 = 0.0 
X3 = -44.5300 
T = 0,0 

SPECIAL FUNCTIONS: 

Y = 2030 . 26*X1_10 . 09*X3 

DERIVATIVES: 

D(X1 /D(T ) = = 

X2 

D(X2 /D(T ) = = 

-2182.68*X1-10.76*X3 
D(X3 /D(T 0 = = 

-12504. 36*X1-62,0*X3 



OUTPUTS : 

TITLE: SES 

TABULATE: T XI X2 X3 Y 

AT INTERVAL 0.300000D-01 
PLOT: X2 

AGAINST: T AT INTERVAL O.lOOOOOD-01 

END CALCULATION WHEN T ,GE. 0.600000 
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APPENDIX D 



The char'a.cteristic equation for the pressure lift only 
model under control conditions was derived from the system 
matrix with the feedback control gains added to the state 
Xg « The result was compared with the characteristic equation 
of Section IV.B.l. 
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The system matrix is 





0 


1 


0 ^ 


A = 


^21 


0 


^23 




^31 


0 


^33 



and with feedback gains j 





0 


1 


0 


A = 


^21 


0 


00 

CM 




^3l"Sll 


+ gl2 


^33‘^^13 



The characteristic equation with gains becomes 

s^-(a33+g^3)s^-(a23g^2‘'^21^® 

■^23^3l“^23®ll 

'^^ 21 ^ 33 ‘^^ 21^13 

It should be noted that the above equation is identical to 
that derived in Ref, 1 except that the gains are shown in 
(4.B.4) 

= -2267.782 
a23 - —11.248 
a33_ = -13002.8 

^33 = 
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From the Riccati equation, it was seen that the gains had the 
following sign values ; 



"■®11 ’ ■^^12 » ''§13 

Substitution of the appropriate signs yields 

s^+(a33+gi3)s2+(a23gi2+a2i)s 

'^^21^33'^^21®13 

"^23^3l"^23^11 

Substituting in the system value parameters, the character- 
istic equation becomes: 

s^+(64,49 3+g3^3)s^+(2 26 7.7 82 + 11.248g3^2)s-11.24 8g^^ 

+ 2267 ,782g3_3 + .5710 
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APPENDIX E 



FIXED AIR FLOW WE STUDY 

A brief study was conducted to examine the effects of 

6 

on the overall frequency response and compared with the op- 
timal control design results. Earlier work indicated that 
the system response was sensitive air flow rate. Therefore 
three values of kg were tested on the pressure lift only 
non linear system with the resultant response curves plotted. 

Values for k_ are : 
b 

kg = 6,144 , Fig. El 

kg = 12.0 , Fig. E2 

kg = 19,0 m Fig, E3 

Note that the high acceleration midband frequencies are 
attenuated as kg increases. The response curve flattens 
out and midband frequencies are reduced by 20%. This result 
would tend to confirm the contention, of Ref. 3, that fil- 
tered fixed air flow rates method of control may have 
advantages. The time response for sinusoidal input remains 
linear throughout the range of interest, The time versus 
acceleration plot shown is a time history for an air flow 
rate of = 60 at an encounter frequency of o)^=45 rad/sec, 

Appendix C presents the DSL Bode program listing. Finally, 
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the system peak frequency magnitude plot for the conditions 
of no control, state variable feedback control and fixed air 
flow rate control is shown in Fig, E8 for Qjq=60. 
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,- 24.00 - 16.00 - 8.00 0. 00 8.00 16.00 24.00 32.00 40.00 48.00 



CRB SES K6=6. 14422 
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o '. 40 
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l '. 60 
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